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Observing  the  Sun  Through  the  Z-mode 


Abstract 

The  directivity  offered  by  the  Z-mode  can  be  very  advantageous 
in  making  radio  astronomical  observations  of  high  resolution. 

In  this  paper  we  have  calculated  the  positions  of  a  1000  km 
circular  orbit  satellite  in  terms  of  the  universal  time  and  the  day  of  the 
year,  from  which  it  can  see  radiation  from  the  Sun  propagating  through 
the  Z-mode. 

*  *  * 

A  satellite  orbiting  in  a  region  where  the  plasma  frequency, 

I  e^N 

l'N  =  y  ,  is  somewhat  higher  than  the  operational  frequency,  f,  can 

detect  radio  waves  from  outer  space  only  through  the  Z-mode. 

The  Z-mode  can  propagate  only  in  the  direction  parallel  to  the 
magnetic  field.  The  direction,  0q  ,  toward  which  we  are  observing 
after  refraction  is  related  to  the  direction  of  the  magnetic  field,  a,  by 


sin  6 

o 


sin  a 


(1) 


eH 

where  fpj  =  is  the  cyclotron  frequency. 

Z  Z 

The  Z-mode  penetrates  the  layer  X  =  1  (i.e.  f^  =  f  )  and  is 

reflected  at  the  layer  X  =  (1  +  Y)  cos^  0q  (i.  e.  fN^  =  (f2  +  ff^)  cos20q)  . 

Z  Z  Z  Z 

If,  therefore,  a  satellite  is  in  the  region  f  <  f.T  <  (f  +  ffu)cos  0  ,  it 

IN  ri  O 

will  receive  only  the  Z-mode. 

Fig.  1  gives  a  schematic  diagram  of  the  Z-mode  propagation. 

For  a  satellite  orbiting  at  1000  km,  a  suitable  frequency  for  observa¬ 
tions  utilizing  the  Z-mode  is  f  01  Z .  Z  Me.  The  ordinary  wave  will  be 

4  3 

reflected  where  the  plasma  frequency  f^  =  Z.Z  Me,  i.  e.  N  =  6  x  10  el/cm  . 
Electromagnetic  radiation  at  Z.Z  Me  propagating  in  the  Z-mode  will  be 


reflected  at  : 
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Electron 
Density 
Profile 


Reflection  Layers 

Extraordinary 

Ordinary 

Z-mode 


Figure  1 


fN2  =[(2.2?  +  fH  2.2]cos20o  .  (2) 

For  mid  latitudes  at  a  height  of  1000  km,  f^~0.9  Me  and  0q  a  14* 

4  3 

which  introduced  in  Eq.  (2)  gives  f^  =  2.  53  Me,  i.  e.  N  =  8  x  10  el/cm  . 

Assuming  a  critical  frequency  for  the  ionosphere  f  »  5.  7  Me, 

5  /  3 

i.  e.  N  =  4  x  10  ei/cm  ,  we  have: 
max  ' 


N 

-  =0.15  and 

max 


N, 


*r 


max 


0.  20. 


From  curves  giving  N/N  vs.  h-h  we  see  that  the  ordinary  wave 

and  the  Z-mode  wave  will  be  reflected,  respectively,  a  few  hundred 

kilometers  above  and  below  the  1000  km  altitude. 

The  earth's  magnetic  field  is  taken  to  be  that  of  a  dipole.  A 

cyclotron  frequency  fH  =  2-8  “7”  =  0.  571  Me  will  be  assumed  at  the  equator 

r 

and  for  the  altitude  of  1000  km. 

2  l/Z 

In  Eq.  (1)  both  f^  =  0.  571  (1  +  sin  $  )  and 
cos  <t> 

(1+3  sin^$)*/2 


sin  a  = 


are  functions  of  the  geomagnetic  latitude. 


'1 
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The  angle  0q,  also  a  function  of  geomagnetic  latitude,  is  obtained  from 
Eq.  (1)  which  takes  the  form: 


sin  6 


\  Ip.  571  (1  + 
y  0. 571  (1  + 


,  •  2,.*l/2 

3  sin  4>)  ' 


3  sin^) 


T727 


2.2  (1  +  3  sin» 


cos  ± 

— rnTr 


<3) 


After  some  trigonometric  and  algebraic  manipulations.  Eq.  (3) 


yields: 


tan  ($-6q)  =  tan  $ 


{- 


1 


.  2,  f\  .  n  (1  +  3  sinz4>)1/z 

sin  $1  +  2  1  +  0.  963  - - 

L  sin  $ 


« 


For  negative  values  of  $,  6q  changes  sign  so  that($-0  )  is 
symmetric  about  $  =  0. 

In  Fig.  2  we  have  plotted  a,  0q,  <4>-©c>) »  and  j9O°-(<j>-0o)J  vs.  $. 
Some  characteristic  values  are  shown  in  Table  I. 


Table  I 


4> 

90* 

O' 

o 

• 

1 

45* 

30* 

0* 

• 

o 

CO 

1 

-45* 

-60* 

• 

o 

O' 

t 

a 

0* 

16.  1* 

26.6* 

41.0* 

90.  0* 

139.0* 

153.4* 

163.9* 

180* 

0 

o 

0* 

9.0* 

14.  0* 

19.4* 

±27. 0* 

-19.4* 

-14.0* 

-9.0* 

0* 

M>-V 

90.  0* 

51.0* 

31.0* 

10.6* 

±27.0* 

-10.6* 

• 

o 

CO 

l 

-51.0* 

-90.0 

0* 

39.  0* 

59.  0* 

79.4* 

117. 0* 
63.0* 

100.6* 

121.0* 

141. 0* 

180* 

Fig.  3  shows  the  angles  $,  a,  0q,  ($-0q)  and  |j)O-($-0o)J  for 
different  cases.  The  angle  ^9O-(^-0q)^J  is  the  angle  that  the  incident 
ray  makes  with  the  earth's  magnetic  axis  directed  toward  the  north 
magnetic  pole. 

As  we  approach  the  magnetic  equator,  the*  penetration  of  the 
Z-mode  will  become  more  shallow,  due  to  its  oblique  indicence.  This 
will  bring  the  reflection  point  of  the  ordinary  and  the  Z-mode  closer 
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Figure  3 


together  and  will  reduce  the  reliability  of  observations  utilizing  the 
Z-mode. 


One  can  undoubtedly  trust  the  Z-mode  for 


which,  as  seen 


from  Fig.  Z  requires  for  the  particular  case  considered  here. 

When  observing  the  sun,  the  angle  =  [9O-($-0q)J  is  the  angle 
that  the  sun's  rays  make  with  the  earth's  magnetic  axis,  measured  on 


the  earth's  magnetic  meridian  whose  plane  contains  the  sun. 


Let  fa)  be  the  angle  that  the  earth  has  advanced  (^l'/day)  on  the 
ecliptic,  starting  from  the  summer  solstice,  and  6  the  longitude  angle 


6 


of  the  midday  meridian  with  the  *<110*E  meridian,  i.  e.  the  meridian 
passing  through  the  magnetic  poles.  The  angle  6  (in  degrees)  is  related 
to  universal  time  (in  hours)  through  Eq.  (5) 

6  =  [l5  x  (Univ.  Time)  -  70*J  .  (5) 

Let  us  also  calltjr  and  the  angles  that  the  sun's  rays  make  with  the 
earth's  rotational  and  magnetic  axes  respectively. 

The  earth's  rotational  axis  makes  an  angle  23*27'  with  the  normal 
to  the  plane  of  the  ecliptic  and  an  angle  11*25'  with  the  earth's  magnetic 
axis. 

The  different  angles  are  shown  in  Fig.  4. 


Figure  4 


From  spherical  trigonometry  we  see: 


cosy  =  cosu)  sin  (23*27*) 

and 

cos^(  =  cost^  cos  (11*25*)  -  sin>|r  sin  (11*25')  cos  6 

Solving  Eqs.  (6)  and  (7)  we  obtain  curves  of  vs.  6  for  different 
values  of  CJ  .  The  results  are  shown  in  Fig.  5. 

The  values  of^(  we  obtain  are  symmetric  in  u>  and  6,  i.e.  we 
have  the  same  values  for  Co  and  (360  -U)  )  and  also  for  6  and  (360  -  6). 
Typical  values  of  are  given  in  Table  II. 


Table  II 


\  u 

6  \ 

0,  360 

30,  330 

60,  300 

90,  270 

120,  240 

150,210 

180 

0 

78.0* 

81.2* 

90* 

101. 4* 

112. 9* 

121.6* 

124.9 

30 

76.  5* 

79.8* 

00 

OD 

LP 

• 

99.  9* 

111.3* 

119.9* 

123.2 

60 

72.  5* 

75.8* 

84.  4* 

98.  0* 

107. 0* 

115. 5* 

118.7 

90 

67* 

70.  2* 

00 

• 

90* 

101.  3* 

109.8* 

113.0 

120 

61.  3* 

64.  5* 

73.0* 

82.  0* 

95.6* 

104.2* 

107.  5 

150 

56.8* 

60.  1* 

68.7* 

80.  1* 

91.  5* 

100.2* 

103.  5 

180 

55.  1* 

58. 4* 

67.  1* 

78.  6* 

90* 

98.8* 

102.0 

In  Fig.  5  we  have  also  included  the  curves  Jj?0- («^»- ©0)]J  vs.  4>  from 
Fig.  2.  Since  both^(  and  [9O-(4>-0oj]  are  the  angles  the  sun's  rays 
make  with  the  earth's  magnetic  axis  we  have  the  horizontal  axis  to  be 
jc  -  po-w-o  q  and  the  two  vertical  axes  are  6  and  <|>. 

The  process  for  finding  the  angle  <|>  that  corresponds  to  a  certain 
6  for  a  given  co  is  shown  for  the  case: 


9 


’A 


6  =  30*  (or  330*)  ,  u>  =  30*  (or  330*) 

where  we  read  4>  =  29*45'  and  4>  =  -14*30'. 

Proceeding  as  above  from  a  larger  scale  diagram,  we  have 
obtained  the  curves  of  $  vs.  6  which  are  shown  in  Fig.  6.  Some 
characteristic  values  of  $  are  given  in  Table  III. 

Table  III 


0,  360 

30,  330 

60,  300 

90.  270 

120,  240 

150,  210 

180 

0.  360 

+3  IjO* 
-13.3* 

+285* 

-15.7* 

+22. 3* 
-22. 3* 

+13. 5* 
-30. 7* 

+  4.3* 
-38.8* 

-45.  5* 

-47.7 

30,  330 

+32.0* 

-12.0* 

+29.7* 
-14. 5* 

+23.  3* 
-21.2* 

+14. 5* 
-29. 5* 

+  5.6* 
-37.7* 

-44.2* 

-46.  5 

60,  300 

+  35.0* 

-  9.0* 

+32. 5* 
-11.  5* 

+26. 3* 
-18.0* 

+17.8* 
-26. 5* 

+  9.3* 
-34. 5* 

+  1.8* 
-41.0* 

-43.  3 

90,  270 

+39.0* 

-  4.  3* 

+36. 5* 

-  7.0* 

+30. 5* 
-13.7* 

+22.2* 

-22.2* 

+  13.7* 
-30. 5* 

+  7.0* 
-36.5* 

+  4.3 
-39.0 

120,  240 

+43. 3* 

+41.0* 

+  1.8* 

+34. 5* 

-  9.  3* 

+26. 5* 
-18.0* 

+18.0* 

-26.3* 

+  11.  5* 
-32. 5* 

+  9.0 
-35.0 

150,  210 

+46. 5* 

+44. 2* 

+37. 7* 

-  5.8* 

+29.4* 

-14.8* 

+21.2* 
-23. 3* 

+14. 5* 
-29.7* 

+  12.0 
-32.0 

180 

+47.7* 

+45.  5* 

-*  38.  8* 

-  4.  3* 

+30. 5* 
-13.8* 

+22. 3* 
-22. 3* 

+  15.7* 
-28. 5* 

+  13.3 
-31.0 

This  formally  completes  the  solution  of  the  problem  but  the 
coordinate  system  in  which  the  solution  is  given  is  not  very  practical. 

To  locate  the  Z-mode  satellite  position  on  a  given  day  of  the  year 
(specified  by  «  ),  at  a  given  universal  time,  we  must  first  find  the 
parallel  circle  where  the  sun  is  located  (the  sun's  latitude,  9,  is  given 
by  sin  9  =  sin  23*27'  cos  w  ).  On  this  circle,  we  then  locate  the  sun's 
position  for  the  given  U.T.  (e.g.  U.T.  =  0,  long.  =  180*). 


n 


On  the  magnetic  meridian  that  passes  through  the  sun's  location, 
take  the  angle  obtained  from  the  curves  of  Fig.  6  where  6*  =  (15*  x  (UT)-70*) 
with  U.T.  given  in  hours.  In  other  words,  the  solution  gives  the  geomag¬ 
netic  latitude  on  the  magnetic  meridian  whose  plane  contains  the  sun. 

For  easier  interpretation,  we  will  convert  this  hybrid  coordinate 
system  to  the  conventional  polar  coordinate  system  of  the  earth. 

Converting  one  system  to  the  other  is  straight  forward  but  it 
requires  some  algebra.  The  two  coordinate  systems  for  a  given  point, 

C,  are  shown  in  Fig.  7. 


P  is  the  North  Pole, 

O  is  the  Center  of  the  Earth, 
A  is  the  North  Magnetic  Pole, 


S  is  the  location  of  the  sun. 

j(os=  X  =  [90-«>-eoO 

cos=  e 

O 


=  4*  (cos  4*  =  sin  23*27'  cos  w  ) 
POA  =  11*25' 

ifpfc  =  6 

SPX  =  180-6  (sin  (S$A)  =  sin  6) 


BOC=  4> 
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From  the  spherical  triangle  SPA  we  have 

sin6  _  sin  PSA 
slnjt  "  sin  1 P25* 


f 


(8) 


and  from  the  spherical  triangle  PSC  (P®A=  PSC)  we  have: 

cos  (F&C)  =  sin  (COO)  =  cos  <P$S)  cos  (c3s)  +  sin  (P$S)  cos  (d3s)  cos  (MA) 
and  therefore 


sin  4  =  cos  4  cos  0Q  +  sin  4*  sin  0q  cos  (PSA.) 


(9) 


Thus  for  given  «  we  obtain  4*  (cos  4*  =  sin  23*27'  cos  <•>  )  and  for 

given  Univ.  Time,  we  obtain  6  (6*  =  15  x  U.T.  -70*). 

Having  w  and  6  we  can  find  y  from  Fig.  5,  4>  from  Fig.  6,  and 
using  the  value  of  $  we  obtain  0q  from  Fig.  2. 

4>  and  are  used  in  Eq.  (8)  to  get  the  angle  PSA,  which, 

along  with  $  and  8q  ,  is  used  in  Eq.  (9)  to  give  us  4  >  the  geographic 

latitude  of  the  point  considered. 

Finally,  from  the  spherical  triangle  SPC  we  have: 


tan 


sin  | 

(p6s  +  cos  -  poc) 

1  sin  1 

(poc  +  cos  -  posY^ 

[ - 2 - J 

l - 2 - ll 

sin 

/  A  A  A  % 

IPOS  +  COS  +  POCJ 

1  sin 

'  p3c  -  c3s  +  p6sU 

1/2 


Calling  "7)  the  longitude  of  the  point  C,  we  see  that: 

A  A  A  A  A  A  /\/\ 

SPD  SSPC  =  EPC  -  EPS  =  EPC  -  6  =  FPD  -  6  =  LPD  -  LPF  -  6=  Tj-110-6 
and  therefore 

a  _ 

SPD  =  ^  -  110*  -  6 
A  /\  A 

Utilizing  the  relations  POS  =  ^ ,  COS  =  0Q  ,  and  POC  =  (90  -  4  )  we  have 


1 
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For  the  given  w  and  U.T.  ,  6,  44  0Qi  and  £  ,  have  been  found 

and  Eq.  (10)  gives  us  ^  which  completes  the  transformation  from  the 
hybrid  coordinate  system  of  6  and  $  to  the  conventional  geographical 
system  of  £  (latitude)  and  (longitude). 

Table  IV  gives  the  values  of  £  and  ^  for  two-hour  intervals 
in  U.T.  and  in  one-month  intervals  starting  (u  =  0)  from  the  summer 
solstice  (  ***  22  June). 

In  general,  there  are  two  possible  satellite  positions  for  solar 
observations;  one  position,  however,  is  usually  quite  close  to  the 
magnetic  meridian,  making  Z-mode  observations  not  too  reliable. 

Thus,  in  Table  IV,  only  the  more  reliable  of  the  two  is  included. 

Fig.  8  gives  the  curves  of  £  for  different  values  of  w  and 
the  positions  of  the  satellite  for  different  U.  T, 

Continuous  lines  give  the  more  reliable  satellite  latitudes  and 
the  dashed  lines  give  the  latitudes  where  the  Z-mode  observations  are 
rather  uncertain.  Magnetic  coordinates  are  shown  in  dashed-dot  lines. 

Because  of  the  many  curve  readings  involved  in  the  calculations, 
which  were  only  carried  to  three  significant  figures,  it  is  believed 
that  the  results  have  an  accuracy  of  about  one  degree. 

Further  accuracy  would  be  of  doubtful  use  since,  by  using  the 
simple  dipole  approximation  for  the  earth's  magnetic  field,  we  have 
introduced  in  our  calculations  an  uncertainty  of  the  same  order  of 
magnitude. 
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Table  IV 


X 

(h) 

0,  360* 

30,  330* 

60,  300* 

90,  210* 

120,  240* 

150,  210* 

180* 

0 

40.  2N 
175W 

37. 7N 

17  5W 

31. 3N 

17  5W 

20. 6S 
176E 

29. 2S 

176E 

35.  7S 

176E 

38.  2S 
176E 

2 

41. 5N 
153E 

39.  2N 

1 53E 

32.  6N 

1 53E 

19.  7S 
148E 

28.  3S 

148E 

34.  6S 

148E 

37.  OS 
148E 

4 

42.  3N 

12  1  E 

40.  ON 

12 1 E 

33.  3N 

12  IE 

19. 2S 
120E 

27.  6S 

120E 

34.  2S 

120E 

36.  6S 
120E 

6 

42.  ON 
88E 

39.  8N 

88E 

33.  IN 
88E 

19.  4S 

91E 

27.  8S 

91E 

34.  3S 

91E 

36.  7S 
91E 

8 

41. ON 
S6E 

38.  8N 

56E 

32.  2N 

56E 

19.  9S 

63E 

28.  6S 

63E 

34.  9S 

63E 

37.  5S 
63E 

10 

39.  5N 
24E 

37.  3N 

24E 

30.  7N 
24E 

21.  IS 

35E 

29.  6S 

35E 

36.  2S 

35E 

38.  6S 
35E 

12 

38. 2N 

4W 

35.  7N 

4W 

29.  2N 

4W 

20. 6N 

4W 

31.  3S 

5E 

37.  7S 

5E 

40.  2S 
5E 

14 

37.  ON 
32W 

34.  6N 

32W 

28.  3N 
32W 

19.  7N 
32W 

32.  6S 

27W 

39.  2S 

27W 

41.  5S 
27W 

16 

36.  6N 
60W 

34.  2N 

60W 

27.  6N 
60W 

19.  2N 
60W 

33.  3S 

59W 

40.  OS 

59W 

42.  3S 
59W 

18 

36.  7N 
89W 

34.  3N 

89W 

27.  8N 
89W 

19.4N 

89W 

33. IS 

92W 

39.  8S 

92W 

42.  OS 
92W 

20 

37.  5N 
117W 

34.  9N 
117W 

28. 6N 
117W 

19.  9N 
117W 

32.  2S 

124W 

38.  8S 

124W 

41.  OS 
124W 

22 

38.  6N 

14  5W 

36.  2N 
145W 

29.  6N 
145W 

21.  IN 
145W 

30.  7S 

156W 

37.  3S 

1 56W 

39.  5S 
156W 

This  accuracy  is  considered  sufficient  as  the  purpose  of  this 
investigation  was  to  locate  the  approximate  positions,  for  a  given  time 
and  date,  from  which  a  2.  2  Me  satellite  radio  telescope  orbiting  at 
1000  km  could  see  the  sun  through  the  Z-mode. 

In  practice,  this  will  allow  us  to  examine  more  carefully  the 
neighborhood  of  these  positions  on  our  satellite  records. 


If  some  increase  in  radiation  is  detected,  more  accurate 
calculations  can  then  be  performed  for  the  exact  time  and  position  of 
the  satellite.  The  electron  density  and  the  magnetic  field,  which  are 
also  available  from  special  instruments  included  in  the  payload  of  the 
satellite,  will  enable  us  to  confirm  proper  conditions  for  Z-mode 
observations.  These  observations  will  give  us  the  means  to  separate 
the  sun's  radiations  from  that  of  other  sources. 

Similar  calculations  can  be  performed  for  Jupiter,  Cassiopeia  A, 
and  other  discrete  sources  of  radiation,  hopefully  allowing  their 
independent  observation  at  long  wavelengths. 
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